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1. Introduction 

For about two decades, nano-structures, like quantum dots, quantum wires or quantum 
wells, are produced by diverse techniques such as etching, local inter-diffusion, particle 
suspension in dielectric media, or by self-assembly in matrices of a host material. They 
display many effects of standard atomic physics by restricting to a confined region of 
space the motion of one to a hundred embedded elementary charge carriers, which may 
be conduction band electrons, valence band holes, or excitons of the semiconducting 
host substrate. In contrast to atoms, two Quantum Dots (QDs) are never identical 
because phonons, surface effects and bulk disorder play a crucial role on their electronic 
properties. But, a QD may be considered as a giant artificial atom, which possesses 
an adjustable quantized energy spectrum, controlled by its size. Therefore, it enjoys 
prospects for an increasing range of future applications: e.g. as a semiconductor laser [1] 
or as single-photon sources [2], as qubits for quantum information processing [3], as 
single-electron transistors in micro-electronics [4], as artificial fluorophores for intra- 
operative detection of tumors, biological imaging or cell studies [5]. 

Thanks to the progress of semiconductor growth technology during the early 
eighties, quantum size effects (QSE) showed up through optical properties of spherical 
semiconductor micro-crystals embedded in an insulating matrix [6,7]. The characteristic 
blue-shift observed in optical spectra of such strongly quantum-confined systems emerges 
in a widening of semiconductor optical band gap, caused by the increasing confinement 
energy for decreasing QD size [8]. It has been also observed in a large variety of other 
confined micro-structures, like quantum ribbons or quantum disks [9], quantum wires 
[10] and quantum wells [11]. A review of empirical and theoretical results on quantum 
confinement effects in low-dimensional semiconductor structures is given in [12]. The 
first theoretical attempt to describe semiconductor QDs has been elaborated upon a 
particle-in- a- sphere model, in the effective-mass approximation (EMA), which assumes 
parabolic valence and conduction bands [8,13-16]. Both electron and hole behave as 
free particles, trapped in a spherical infinite potential well, and move with their common 
effective masses in the considered semiconductor. The electron-hole Coulomb interaction 
is included, and the excitonic contribution to the QD ground state energy is taken into 
account by Ritz' variational principle. Other EMA models have been built upon finite 
potential wells, improving agreement with experimental data for a significant range 
of QD sizes [17-21]. In addition to spherical clusters, the case of cylindrical shaped 
micro-crystals has been carefully studied [17,22-24], as well as the case of quantum 
wires [25,26]. More sophisticated models, which consider non-parabolic valence and/or 
conduction band(s), have been also developed [27-31]. Modern approaches and how 
they can be applied to large structures and compared to experiments are discussed in 
the book [32] as well as in the recent review article [33]. 

Among many important topics, it is the physics of atom-like behavior of QDs, which 
is nowadays most vigorously investigated in quantum confined systems, for potential 
technological applications. Of particular interest is the interaction with an ambient 
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electromagnetic field, giving rise to the so-called quantum-confinement Stark effect 
(QCSE), which has been studied for example in CGLSo.i25eo.88 [34], in CdSe [35] or in 
InGaAs micro-crystals [36]. It consists of an observable red-shift of an optical transition 
induced by the presence of a constant external electric field [37-40]. In recent years, 
some works have also dealt with ac-electric field [41,42]. Stark effect leads to an energy 
shift of the exciton photoluminescence as well as a corresponding enhancement of its 
recombination lifetime [43]. The electric field dependence of QCSE was first studied 
in GaAs — AlGaAs multi-layers quantum wells [44]. Exciton energy shift peaks were 
experimentally observed and successfully compared to theoretical results obtained by a 
perturbation method introduced in [45,46], as the applied electric field is perpendicular 
to the plane of the layer wells, within which the electron-hole Coulomb interaction is 
negligible. However, in spherical QDs, the Coulomb potential turns out to be more 
important, and cannot be discarded [34,47]. 

Over the years, the spherical shape of QD has remained popular in the study of 
QCSE [48-52]. But, to the best of our knowledge, no simple comprehensive model, 
which describes Stark effects in spherical semiconductor nano-structures with analytic 
results, has been worked out. In this paper, we propose to use the EMA model for 
spherical micro-crystals, to establish analytically some criterions on the QD radius and 
on the electric field amplitude, and to understand why presently known results fail so far 
to correctly describe QCSE for a wide range of QD radius. To this end, in Section [2} we 
shall introduce the electric field free model first, and recall some of its general properties. 
The next two Sections [3] and [4] are devoted to the analysis of Stark effects in spherical 
semiconductor QDs, first with the inclusion of electron-hole Coulomb potential and 
second with an additional polarization energy. In the concluding section, we summarize 
our main results and indicate possible future research perspectives. 

2. EMA Quantum Dot model 

A standard EMA model with infinite spherical confining potential well, without electron- 
hole spin coupling and external magnetic field, allows to perform analytically most of 
the calculations on spherical semiconductor nano-structures interacting with a fixed 
external electric field. There exists other models with parabolic confinement [53, 54] 
or parabolic potential superimposed to an infinite potential well [55], but the concept 
of a QD size is then not so well defined. As Stark effect in semiconducting micro- 
crystallites manifests itself through an energy shift of the electron-hole total energy 
levels, we have to deal mainly with energy eigenvalue differences of a Hamiltonian. 
However, two different energies have to be computed within the same theoretical QD 
model. Even if this model does not fully describe the QD behavior in the absence of 
electric field, particularly for small QD radii, it can be still used, since it gives rather 
satisfactory theoretical predictions on Stark effect. The overestimation made for the 
electric field free electron-hole pair energy levels should also appear in the interaction 
of the electric field with electron-hole pairs [20]. Then, despite intrinsic limitations, an 
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approximate value of the Stark shift can be obtained under some consistency conditions. 

2.1. Consistency conditions 

Since most synthesized nano-crystallites possess an aspect ratio (defined as the ratio 
between the longest and shortest axes of the QD) smaller than 1.1, even if higher 
aspect ratio micro-crystals would exist, the hypothesis of a QD with spherical symmetry 
appears often as quite reasonable. 

Realistically, the effective potential at the QD surface is finite, and has a standard 
amplitude from 1 to 3eV [20]. This value justifies already the use of an infinite confining 
potential well, as mentioned earlier, because it is generally quite large as compared to 
typical electron and hole confinement energies usually involved, which increase as cx R~ 2 
for decreasing QD radius R [14]. Therefore, the tunnel conductivity through the QD 
boundary is vanishingly small, except of course for very small QDs. Futhermore, the 
infinite potential well approximation implies that charge carriers inside the cluster are 
insensitive to its outside surroundings, particularly to any externally applied field, as 
far as considerations on QCSE are concerned. Although the surrounding effects may be 
sufficiently small to be neglected, the presence of a large external field can significantly 
modify the inside behavior of the micro-crystallites. Thus, the electric field amplitude 
outside the QD should not then exceed a threshold, fixed ad hoc by the height of the real 
confining potential step. This constraint should be referred to as the usual weak electric 
field limit. An inequality, which analytically expresses its validity by linking the electric 
field amplitude to other physical parameters of the problem, is to be determined later. 
It allows to evaluate an approximate value for the maximal electric field amplitude to 
apply on the QD, while respecting the weak field limit. 

Lastly, for small nano-crystals of typical sizes of less than a hundred lattice spacings, 
there exist magic numbers for which clusters remain stable: e.g. crystalline silicon only 
stay coherent as clusters of 5^12, Si 33 , Si 3Q and Si^, if they contain less than 60 silicon 
atoms [56]. Their band structure are so deformed such that it becomes impossible to 
use the parabolic shape of conduction and valence bands, required by EMA models. 
However, if there is no potential well — i.e. if there is no semiconductor micro-crystal 
embedded in the surrounding insulating matrix — , no electrons should be excited, and 
no holes should appear. Stark effects must then vanish in small QDs for any electric field. 
In the weak field limit, it has been shown that, in semiconducting rectangular quantum 
boxes, the Stark shift of the confined exciton ground state presents three contributions, 
each of them going as the fourth power of an edge length [40]. Hence, for a spherical 
potential well of radius R, the Stark shift is expected to scale as cx i? 4 . 

2.2. General considerations 



written in spherical coordinates (r,9,<p). According to the EMA model, the total 



We work with the infinite confining potential well V(r) = V(r) 
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electron-hole pair Hamiltonian operator H reads, in units where H = 1 
H = H C + H h + V c (r ch ) = -^5- - + V(r e ) + V(r h ) 



2m* 2m£ ftr eh 

where /t = 47re, e denotes the semiconductor dielectric constant, r e h = |r e h| = |r c — r^l 
the electron-hole relative distance, m* h the effective mass and if e ,h the confinement 
Hamiltonian respectively of the electron and of the hole and Vc(r e h) the electron- hole 
Coulomb interaction. Without loss of generality, the semiconductor energy band gap 
Eg is set equal to be zero for convenience. 

In the absence of Coulomb potential, electron and hole are treated as decoupled 
particles, the QD wave function should be then factorized into separable electronic and 
hole parts \l/(r e ,rh) = ^(r^^^h). The orthonormal eigenf unctions ipi nm ( r ) are labeled 
by three quantum numbers ZgN, nGN* = N \ {0} and me [— I, /]. If < r < R, 

where Y[™(#, </?) is the spherical harmonic of orbital quantum number I and azimuthal 
quantum number m, ji(x) the spherical Bessel function of the first kind of index I and 
of variable x, and {ki n } ln the wave numbers set, defined as the n th non-zero root of 
the spherical Bessel function }i(x) thanks to the continuity condition at r = R — the 
presence of an infinite potential imposes that ipinm( r ) = 0, if r > R [13]. The respective 
energy eigenvalues for electron and hole are expressed in terms of {ki n }i n as 

k 2 

111 ~ 2ml h R 2 ' 

The continuum density of states of the semiconductor bulk should show atomic-like 

discrete spectrum with increasing energy separation as the radius decreases. 

Because of the explicit micro-crystallites spherical symmetry breakdown in the 

presence of a Coulomb potential, the exact determination of eigenfunctions and energy 

eigenvalues for Eq. ([T|) is arduous. Treating the interplay of the Coulomb interaction, 

scaling as oc -R -1 , and the quantum confinement, scaling as oc R~ 2 , constitutes the 

common approach to this problem. To handle these competing contributions, two 

working regimes are singled out, according to the ratio of the QD radius R to the 

Bohr radius of the bulk Mott-Wannier exciton a* = -r— , \i being the reduced mass of 

e 2 jJL 

the exciton [16]. In the strong confinement regime, corresponding to a size R < 2a*, 
the potential well strongly affects the relative electron-hole motion, and exciton states 
consist of uncorrelated electron and hole states. The weak confinement regime, valid 
for a size R > 4a*, leaves the electron- hole relative motion and its binding energy 
unchanged. The exciton character of a quasi-particule of total mass M = m* + m£ is 
conserved. As its center-of-mass motion remains confined, it should be quantized [16]. 

Even if we focus on Stark effect in the strong confinement regime, we shall briefly 
present the consequences of the previous single EMA model in both strong and weak 
confinement regimes. Despite its simplicity, this model seems to be able to apprehend 
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correctly the QCSE, at least for a range of sufficiently small QD sizes, and to yield 
numerical values, which agree with experimental results. 

2.3. Considerations on strong and weak confinement regimes 
2.3.1. Strong confinement regime 

In this regime, the Coulomb potential is treated as a perturbation with respect to the 
infinite confinement potential well in a variational procedure, which shall be extended 
to the case of an applied electric field. The ground state energy of the exciton should 
be evaluated with the following trial wave function 

0(r c ,r h ) = ^oio(r c )^oio(rh)0rci(rch), (3) 
with rc i(r c h) = 0rci(^eh) = e~^ rch , where a is the variational parameter, r ejh = |r e h | and 

Sm &e,h) 



^oio(r e ,h)=^oio(^e,h) 



The variational wave function of Eq. ([3]) implies that both electron and hole should 
occupy primarily their respective ground state in the confining infinite potential well, 
as described by the product ^oio( r e)V ; oio( r h)- It should also exhibit, via the function 
re i(r e h) of the relative coordinates r c h, an exciton bound state behavior, analogous to 
the ground state of an hydrogen-like atom with appropriate mass /i and Bohr radius a*, 
up to a normalization factor, especially if <t _1 oc a*. 

Despite the breakdown of translational invariance of the Coulomb interaction by 
the spherical confining potential, Fourier transform formalism in relative electron-hole 
coordinates allows to establish integral representations for quantities such as the square 
of the norm {(j)\4>) of the trial function 0(r e , r^) or the corresponding Coulomb potential 
diagonal matrix element (0|Vc(r e h)|0) 

1 ffdxdy 



ffdxdy g . n2 sin y^ sin U aRx \ e -°Rv 
JJv % y 



\Vc(r A )\<t>) = -~ [[ — ^sin 2 (vrx) sin 2 (7ry) sinh^^e"^, 
kRo-JJ v x y 

where T> = {(x, i/) 6 R 2 / 0<x<y<l}. A Taylor expansion of expressions Eqs. 
with respect to the dimensionless parameter aR near zero yields 

= 1 - BaR + 0(a 2 R 2 ), 
(0|^c(r e h)|0) = ~ {A - aR + 0(a 2 R 2 )}. 



(4) 



Thus, an expression of the mean value of the total Hamiltonian Hq in the strong 

il pa 
1 

Ilia" 



confinement regime in terms of a dimensionless variational parameter a', defined by 
a' = era*, and of the binding exciton Rydberg energy E* = - — — can be obtained as 
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where the correction terms "... " go to zero as soon as — goes to zero, and where 

a* r-i 

i?eh = + -^01 * s the electron-hole pair ground state confinement energy l§|. The 
variational parameter a' is determined to be a' Q = AB' « 0.9956, to minimize the value 
of the electron-hole energy 

,2 



ET ng = E ch -A^--AB' 2 E\ 



e 
~~kR 

This formula has been already analytically obtained with trial functions showing the 
same global form as 0(r e , r h ), but with an interactive part equal to re i(r e h) = (1 — § r e h), 
instead of re i(r e h) [16,57]. It is obvious that re i(r e h) consists of the two first terms 
of the Taylor expansion of </> re i(r e h), in the limit |r c h < aR <C 1. Because of the 
infinite confining potential well assumption, the total excitonic energy is overestimated 
in comparison to experimental data for small QDs. A successful method to subtract off 
this over-estimation consists in adopting a model in which a confining finite potential 
step of experimentally acceptable height is restored [20]. 

2.3.2. Weak confinement regime 

In this regime, electron-hole pair states consist of exciton bound states. The Coulomb 
potential and the kinetic energy in the electron-hole relative coordinates are of the same 
order of magnitude because the QD size allows a partial restoration of the long range 
Coulomb interaction between the charged carriers inside the QD. Then, the essential 
contribution to the ground state energy of the exciton is —E*, the ground state energy 
of a hydrogen-like atom mass /x. Furthermore, the total translational motion of the 
exciton, thought as a quasi-particle of total mass M, should be restored and contributes 

to the exciton total energy by 2j\^^2 ' ^ e typical kinetic energy term of a free particle 
confined in a space region of size R. As a first approximation, the ground state energy of 
the exciton trapped inside the QD is then the sum of these two energetic contributions, 
but this is not totally satisfactory. To improve phenomenologically its accuracy in regard 
to numerical simulations, a monotonic increasing function 77(A) of the effective masses 

772* 

ratio A = — - has been introduced in [16], and has been inserted into E^^ as follows 

pwcak rp* j 

eh 2M(R-r](\)a*) 2 ' 

The QD size renormalization term 7?(A)a* is a dead layer [58]. Although it could be 
successfully described as a quasi-particle, the exciton is not itself an indivisible particle. 
Its center-of-mass, whose motion is quantized, could not reach the infinite potential 
well surface unless the electron-hole relative motion undergoes a strong deformation. 
The picture of a point-like exciton should be dropped in this region of space. The 
exciton should be preferentially thought as a rigid sphere of radius r)(\)a*, where 77(A) 
is numerically determined to get a better fit of experimental results [16]. 

§ I Appendix A| compiles a register of all constants, which appear in this paper. 
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3. QD Stark effect without polarization energy 

The diagonalization of the Hamiltonian of an electron-hole pair trapped in an infinite 

potential well under the influence of an external constant electric field is, in principle, 

an exactly solvable problem. Even if the presence of a constant electric field explicitly 

breaks both spherical QD symmetry and electron-hole Coulomb potential translation 

invariance, in a spherical QD, we shall consider an applied electric field E a along the 

direction z of a three-dimensional cartesian coordinates system with its origin located 

at the QD center. This is not the cases of quantum wells [45,46], quantum wires [59], 

quantum disks [39], rectangular [40], cubical [60] or confined by parabolic potential [61] 

quantum boxes, in which the electric field direction plays a significant role. As the 

inside semiconducting QD dielectric constant e is larger than the outside insulating 

matrix dielectric constant e', the electric field Ed inside the QD is different from E a . It 

E a 

is given by Ed = -, ^ , where g is a geometrical depolarization factor, which 

(i -g) + g£r 

equals - for a sphere, and e r = — is defined as the relative dielectric constant [62]. 
The dielectric constant difference also implies the existence of a polarization energy 
term P(r c , rh), introduced in [8], which shall be neglected in this section, but taken into 
account in the next one. This hypothesis allows to study in more details its relative role 
vs. the Coulomb potential, because they both scale as oc R~ x . 

Let us define the electron and the hole (of respective electric charge =pe) interaction 
Hamiltonian with the electric field Ed, in spherical coordinates, as 

W e ,h(r e ,h) = ±eE d • r 0>h = ±eE d r Cjh cos 6> e , h , (5) 

where E d = |Ed| is the electric field amplitude inside the micro- crystal. As the function 
0(r c , rh) does not provide any further contribution to the excitonic energy in the presence 
of the electric field, i.e. (<f)\W e (r e )\<p) = — (<p\ Wh(i"h)|0), an appropriate form for the trial 
wave function should present some other dependence on the electron and hole space 
coordinates to be determined later. 

3.1. Justification of the variational trial wave function form 

To apprehend the effect of the induced electric polarization, we follow a reasoning 
made in [45], and study the interaction between the charge carriers with the ambient 
electric field, neglecting the Coulomb potential. To this end, we define the individual 
Hamiltonian H' c h of a confined electron or of a confined hole interacting with Ed as 

H' ch = H Bih + W e ,h(r e ,h). 

As mentioned as a consistency condition in Subsection \2.1l we can assume that the 
electric field amplitude is sufficiently small so as to consider the Hamiltonian interactive 
part W e ^(r e ^) as a perturbation to the confined Hamiltonian H e ^. In this weak field 
limit, the absolute typical interaction energy of the electron or of the hole under the 
electric field influence E e \ e = eE d R should be treated as a perturbation compared to 
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their typical confinement energy E ch = Eq\, i.e. E e \ e <C E Cth . In real atoms as well 
as in quantum wells [46], in quantum wires [59] and quantum boxes [60], Stark shifted 
levels show a typical quadratic dependence on the electric field amplitude, a similar 
behavior is expected here. In the following, to justify the form of the variational wave 
function leading to the QCSE, we investigate perturbative results on decoupled confined 
electron and hole, interacting individually with an external electric field, by performing 
first, a standard second-order stationary perturbation theory, and second, a variational 
procedure. 

A second-order perturbation computation of the Stark shift undergone by the 
electron or by the hole ground state, yields 

A£ c S r k pcrt = -r pcrt m c * h e 2 £X + 0(£ d 3 ), 

32 . k 2 
where the constant r Dert = — n 2 } —r. — ln ~ 0.01817. 

3 ^ [ki n — 7T 2 ) 5 
n>l v ln ' 

In order to account for the electric field direction along the z-axis in the variational 
principle, the trial wave function should show a deformation away from the spherical 
shape, which squeezes or stretches the electron or the hole probability density along 
this particular direction. The variational trial function is chosen of the form $ e ,h(r e ,h) = 
^oio(r c ,hVc,h(rc,h), where y?e,h(r c ,h) = e^^ 7 " 6 * 008 " 6 - 11 . The variational parameters a e>ix 
have the dimension of an inverse length so that, in the weak field approximation, we 
can assume that o~ e ^,R <C 1. The difficulty in this problem is in the calculation of 
the square of the norm of the trial function $ e ,h( r e,h)- However, it admits an integral 
representation, on which a Taylor expansion in the neighborhood of the dimensionless 
parameter a e ^R = can be performed 

($ e ,h|$e,h) = „ / — sin 2 (7nc) sinh(a 0i h-Rx) 

= i + f< h ^ + o« h ;? 4 ). 

The mean value of the confinement Hamiltonian if e ,h in the quantum state defined by 
the trial function $ e ,h(i"e,h) is exactly determined as 





He,h 


*e,h> 




$e,h> 



E 



e,h 



and the mean value of the interaction Hamiltonian W e h( r eh) as 











$e,h> 



Then, the total Hamiltonian H' c h mean value, up to the second order in a c ^R, is 

0-2 c 

E e v, H : eE c \R <Tp h + . . . 

* 8m* eh 3 





H 'c,h 
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A choice of the variational parameter a® h = — m* h eEdR 2 gives the ground state energy 

3 

E' c h of the confined charge carriers in interaction with the electric field, and the Stark 
shift by subtracting their respective ground state energy, as follows 

A-k'sttrk var = E' h - E Cjh = -r var m* >h e 2 ^-R 4 , 



2C 2 

where r var = — — ps 0.01776. 

These two Stark shift expressions present the same dependence on physical 
parameters: they both scale as oc m* h e 2 E 2 R A . The Stark shift contribution is clearly a 

second order term in the dimensionless parameter — — <C 1, with respect to the electron 

E e ,h 

or the hole confinement energy i? e ,h- The difference between the previous methods is 
quantifiable by evaluating the relative error between the values of the proportionality 
constants r pcrt and r var , which is ~ 2%. This small relative error supports the validity 
of the new trial wave function $(r c , r^) in the presence of the electric field, defined as 

$(r e , r h ) = 0(r e , r h )( / j c (r c )v?h(r h )- (6) 

Since, the function $(r c , rh) has the part </>(r e , rh), describing the electron-hole Coulomb 
interaction both occupying the ground state of their respective confinement Hamiltonian, 
and the electric field interactive part </9 e ( r e) ( / J h(i'h), liable for the individual electron and 
hole behaviors in the ambient electric field Ej. 

3.2. General results on Stark effect in semiconductor Quantum Dots 

As already mentioned, we add the interaction Hamiltonians We^r^h) between the 
charge carriers and the electric field to the QD model Hamiltonian Hq, introduced 
in Section [21 in order to apprehend QCSE in spherical semiconductor micro-crystals 

H = H + W c (r c ) + W h (r h ). (7) 

In the weak field limit, the variational procedure is to be applied using the trial function 
$(r c ,rh), introduced in Eq. (jBJ). To this end, we use a reasoning similar to that of 
Subsection \2.2\ i.e. the Fourier transform formalism in the relative coordinates can be 
used once again quite advantageously. This formalism leads to integral representation 
of the square of the norm of the trial function $(r e , rh) and of the mean value of the 
Coulomb interaction matrix element in the corresponding quantum state 

1 f 1 ff dx dy 



[ 1 ' R 2 3a 



a 



f ff dx dy 

/ d£ / / sin 2 (7nr) sin 2 (ny) 

J-i JJv x y 



j sinh(p e (£)aRx) e - p ^ )aRy smh(p h (£)aRx) e~ pe ^ R y 

(^c(r ch )|$> = -44; I ^11 --si^Msin 2 ^) 
kR aR J_ ± JJ V x y 

smh(p e (£)aRx) e - ph ® aRv smh(p h (£)aRx) e - p ^ )(TRy 



(8) 
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l_2^£+^,for-l<£<l. 
a a z 

Direct calculations prove that Eqs. dHJ are valid if and only if the variational 
parameters a and cr c h satisfy the inequality 

< e ■ a C)h < a, (9) 

where e = exp(l). Hence, Eq. (Q is a consistency condition, which analytically 
determines the range of acceptable electric field amplitudes. First, following the 
variational results for the interaction between the electric field with the electron or with 
the hole, we expect that cr c h oc m* h e£'di? 2 . Second, following the variational results 
on the electric field free interactive electron-hole pair through the Coulomb potential, 
we also expect that a oc a* -1 . Then, after trivial rearrangements, we remark that 

— — oc — , i.e. the weak field limit should remain valid if — — does not exceed, up 
E e , h a a* E e>h 

R 

to a dimensionless proportionality factor to be given later, the order of magnitude of — , 

a* 

the ratio which characterizes the strong confinement regime. The charge carrier energy, 

in the presence of the electric field, should be at most of the same order of magnitude 

of a first term correction to their confinement energy in the strong confinement regime, 

which corresponds to the absolute value of the typical electron-hole Coulomb interaction 

energy, because they both scale as oc Rr 1 . 

In the limit of vanishing electric field, i.e. in the limit — > 0, the expressions in 

a 

Eqs. (JHJ) allow to retrieve the expressions for the square of the norm and for the Coulomb 

potential mean value without electric field expressed by Eq. (j3j). Moreover, in the weak 

field limit, based on the decoupled electron- hole point of view presented in Subsection 

\3. 11 we expect that the Stark shift for the coupled electron-hole system should scale as 

R 2 

oc (m* + mV)e 2 E\R oc E eh — -. Then, to get at least the lowest order contribution to 

a* 2 

this Stark shift, it is necessary to perform a Taylor expansion of the total Hamiltonian 

H mean value up to the second order in the variational parameters. However, as we 

shall see in the following subsection, this first contribution is not sufficiently accurate to 

fit experimental data, because it does not account for the electron-hole coupling through 

the Coulomb interaction. This is the reason why we shall carry on the expansion up to 

R 

the third order, since we will also get the first correction in — to the Stark shift, which 

a* 

expresses the presence of the Coulomb potential in the strong confinement regime. 

Finally, we obtain the interaction Hamiltonian W e ,h( r e,h) mean value from the 
square of the norm of the wave function $(r c ,rh) by taking its logarithmic derivative 
with respect to the variational parameters o~ c h 

As we have build the trial function <3>(r e , i"h), so that it possesses the properties of 
the functions <f)(r c , r^) and <3> e h(r e h), the exact mean value of the electric field free 
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Hamiltonian H should lead to 

(swig) _ F ,^,^L,^L m) 

($|$) eh+ 8/i + 8m c * 8m*' 1 J 

where the electron-hole pair ground state confinement energy E e ^ is provided by the 

2 2 

function ^oio( r e)V'oio( r h), and the contributions — and — to the total kinetic energy 

8{j, 8m* h 

are respectively due to the Coulomb potential and to the interaction between the charge 
carriers and the electric field. However, in addition to these four expected terms, direct 
calculations exhibit a further contribution to the mean value of H Q of the form 

K(a,a c ,a h ) = - j — + — 1|_ - _ j __ 

Real physical quantities should be invariant under the electron-hole exchange symmetry 
defined by the exchange of their coordinates, their masses and their electric charges 

r c ,h -> r hie , m* h -> m£ c , e -> -e. (12) 

As we expect that cr Cj h oc m* ch eE&R 2 , under an electron-hole exchange, the variational 
parameters should transform as a c ^ — > — a^ c . Therefore, the trial function $(r e , r^), 
the confinement Hamiltonian H , the Coulomb potential Vc(r e h) and the interaction 
Hamiltonian W e (r e ) + Wh(i*h) should remain invariant under the electron-hole exchange, 
as well as the norm of $(r c , r^) and the mean value of these operators. But, the further 
contribution will not, since it changes as K(a, a e , a^) —> —K(a,a e ,a\ l ). Because of the 
mean value of the confinement Hamiltonian H invar iance, it should not bring any new 

($|#o|$) 

contribution to real physical quantities, and should be discarded from — /$|$\ — ' as 
given by Eq. ffTTl) . 

3.3. Stark effect in strong confinement regime 

As there seems to be no way to analytically compute the integrals in Eqs. (jSJ), we propose 

to Taylor expand them in the strong confinement regime, i.e. when aR C 1. To perform 

this expansion, we have to specify the QD radii region, in which the following expressions 

are valid. For this, we shall assume that the quantities p e ^(£)o~R in the arguments of 

the functions exp(x) and sinh(x), appearing in Eq. (jSJ) should be sufficiently small. As 
3 

Pe,h(0 < _ f° r ai l £ ^ [ — 1) 1]) we wm only consider the range of QD radii such that 
2 

R < Tr~, a " \ so that p e h(£)o-Rx, p e h(£)aRy < 1, if — 1 < £ < 1 and < x < y < 1. 

OCT 

Thanks to the consistency condition Eq. (jHJ), the variational parameters a and o" e ,h can 
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be of the same order of magnitude. Then, we deduce, up to third order in o~R, that 

($|$) = 1 - BoR + Ca 2 R 2 - Da 3 R 3 

2 2 

+ {C- D'aR f c+ Q ah R 2 - D"aa c a h R 3 + 0(a 4 R 4 ), 

(mc(r eh )m -e 2 {A + B>aR + CV2jR2 + ^ + al)R2 



+C 2 a c a hJ R 2 + 0(a 3 i? 3 )}, 

(*|w. (r «) + w M _ _ eE ^ + CTh)Jj2 r o _ CVfl + 0(<j2fl2; 



($|$) uv ' {3 

The mean value of the total Hamiltonian H, under the influence of the electric field on 
both electron and hole is then expressed as an expansion in powers of the variational 
parameters, up to third order in aR, 

« = E , h - Ail + + ^L + ^L_ 2BW 

($|$) k.R 4 8m* 8m* 

- 2C'E*—a' 2 - C a ° + a ^ — - C <7G<7h — 

a* 1 /i a* 2 /i a* 

- e£ d (a e + ^h)i? 2 {| - C "'^'}+ • • • ( 14 ) 

We minimize the previous matrix element with respect to a' and a c ^ to obtain an 

approximate value of the ground state total energy. Their values are determined to the 
R 

first order in — to insure the coherence of the expansion we made as 
a* 



R 1 % p 2 F 2 f? 4 f? 

a' Q = AB'{ 1 + 8C'4 - *-CC"{ml + m£)^f 



4(7 

<h=^<h^d^l + 4 



a* ) 3 v c n/ £* a* 



a* 



(15) 



n z n c 

The Stark shift is identified as terms scaling as oc E 2 , andjj] 

A<T = -r var (m: + m*)e 2 £ d 2 i? 4 {l + 81^}, (16) 

where r var appears as of universal character, while the constant depends on the 
semiconductor. In terms of the effective masses of the electron and the hole, it can be 

|| From Eqs. (| 13[) . we confirm that K (<r, a e , Ch), the contribution discarded from the confinement 
Hamiltonian Hq mean value because of the electron-hole exchange symmetry, does not contribute to 
Stark effects 

K(a,a e ,a h ) = 

I ^ ■* J I '"h '"e J 

Such terms contribute to the total Hamiltonian mean value H up to the third order in — . These third 

a* 

order terms just contribute to the variational parameters a' and a® , up to the first order, but not 
at all to the electron-hole pair ground state energy. Therefore, the last possible contribution to Stark 
effect should come from K(a' , a® , o^), where we replace the different variational parameters by their 

4C 

respective zeroth order expressions a' m AB' and ct° h ss -r-wi* ^eE^R 2 , according to Eqs. (Ti5j) . Then, 

1 o 

we straightforwardly verify that K(a' , a®, cr?) = 0, which is what was expected. 









CTe 1 






IK 


m* J 
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expressed as 



ml m* J C 



The first contribution to this shift is simply the sum of the Stark shift contributions 
undergone by the ground states of both electron and hole taken individually as computed 
in Subsection \3.1\ Because of the dependence of the constant on the effective masses 
the second contribution to AEf™^ indicates the existence of a further coupling 



between the electron and the hole, which appears as a standard dipolar interaction 
between two opposite electric charge carriers. This interpretation is a question of point 
view. Until now, we have considered that the interaction between the electron or the 
hole with the external electric field takes place individually, whereas they interact only 
through the Coulomb potential. This physical description justifies a priori the validity 
of the strong confinement regime assumption, for which the exciton states consist of 
uncorrelated individual confined electron and hole states. It then allows to intuitively 
build a coherent model in order to describe the QCSE in QDs in this regime, and also 
simplifies the calculations in practice. In spite of these advantageous properties, the 
previous remark suggests that this picture should be revised. 

Actually, the total Hamiltonian electric field interaction part W e (r e )+Wh(rh) should 
also be written as 

W(r ch ) = W c (r c ) + W h (r h ) = E d • d ch , 

where d c h = er h is the exciton electric dipole moment. This is the standard dipolar 
interaction Hamiltonian of an electric dipole. It satisfies the electron-hole exchange 
symmetry, while the individual interaction Hamiltonian W e ,h(r e ,h) transform themselves 
one into another. In the strong confinement regime, despite the importance of 
confinement effects on excitonic ones, the dipolar interaction point of view expresses the 
remnant of electron-hole pair states, thought as exciton bound states under the influence 
of the electric field. It suggests the inclusion of a further term in the Hamiltonian 
H describing the exciton-electric field interaction, which accounts for the polarization 
energy of the electron-hole pair, due to the difference between the dielectric constants 
of the semiconductor QD and the surrounding insulating matrix. 



3.4- Comparison with experimental data 

In order to test the relevance of our model, we shall compare our predictions to real 
experimental data given in [34] and to other computational data of [47] in spherical 
CdS .i2Se . 88 micro-crystallites with material parameters: e = 9.3, m* = 0.13m e , 
m^ = 0.46m e , E* = 16meV and a* = 49A, where m e is the electron bare mass. The 
electric field amplitude inside the micro-crystal is fixed at E^ = 12.5kV.cm _1 . We note 
that an earlier work [51] has used numerical diagonalization of the total Hamiltonian H, 
with the same material parameters for spherical CcfSo.12Seo.88 QDs as well as electric 
field, to obtain theoretical predictions. 
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Figure 1. Stark shift for confined interactive electron- hole pair as a function of the 
QD radius including the Coulomb interaction and excluding the polarization energy 
up to the zeroth ( — ) or to the first ( — ) order in comparison with experimental results 
(+) [34]. « -0.1629. 
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3.4-1- Case of real experimental data 

Figure [1] presents a comparison between results we obtain and experimental values for 
spherical CdS'o.i2>S'eo.88 micro-crystallites [34]. Two exciton peaks are experimentally 
resolved, which are attributed to the transitions from the highest valence sub-band and 
from the spin-orbit split-off state to the lowest conduction sub-band, with an energy 
splitting is about 0.39eV independently of the QD radius [34]. The experimental values 
depicted by crosses in Figure [T] consist of mean values of the Stark shift of these two 
types of excitons. They seem to indicate that the Coulomb interaction is sufficient 
to explain correctly the amplitude of the Stark effects experimentally observed, as we 
expect, in the range of validity of QD radii. 

In the strong confinement regime, our approach offers a model capable to describe 

R 2 R 
QCSE at least for QD sizes — < — 7 . But, a' n is itself a function of — , which is 

a* 3a a* 

still considered as a small dimensionless parameter in the strong confinement regime. 

Here, the part of a' which depends on the electric field is negligible, because it scales 

as oc m c e 2 E\R A . This is at least of the same order of magnitude as the exciton 

Rydberg energy as soon as R < 50A, if the electric field amplitude is fixed at 

CC"(m* + m*) 

E d = 12.5kV.cm _1 , and - ° — — ~ 0.0552m e . Hence, according to Eqs. (TT5]) . 

1 Ad (_y 

a' m AB' < 1 + 8C' — >. Therefore, up to first order in — , our predictions should be 
{ a* ) a* 

valid for QD radii 

R < — — r w 0.6080a*. (17) 

~ 2(3B> + AC) V ' 

According to this effective constraint, in the case of Cd5o.i2*Se .88 micro- crystals, this 
approach should lead to acceptable results in regard to experimental data as long as 
the cluster radius does not exceed 30A. Figure [T] shows that the absolute value of the 
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Stark shift, computed up to the zeroth order, is significantly overestimated, except for 

a minor range of small QD radii compared with the one we expect. The results become 

much more accurate, if the Stark effects are computed up to the first order. In this case, 

Figured] exhibits a good agreement with the experimental data over the whole expected 

region of micro-crystals radii. In this domain of validity, the first order calculation seems 

to be efficient enough to describe QCSE in spherical semiconductor QDs. As soon as, the 

QD radius exceeds the maximal value for which the strong confinement regime is valid, 

our results diverge significantly from experimental data. Furthermore, Figure 4 in [51] 

shows Stark shifts of the ground state and the two first excited states of the confined 

exciton, obtained by numerical diagonalization of the Hamiltonian H . Agreement with 

experimental data from [34] is also reported in the ground state case. Whereas [51] 

gives a range 20A < R < 100A, in its validity domain R < 30 A, we may say that our 

analytical approach is totally consistent with the numerical approach of [51]. 

To determine the maximal electric field amplitude, for which the weak electric field 

limit assumption remains valid, the consistency condition Eq. ([9]) is reconsidered for 

the electron or for the hole, in which the respective variational parameters are replaced 

by their variational values. After summing the expressions for the electron and for the 

R 

hole, up to first order in — , we deduce that 

a* 

E ] 11 _ 
— — < ps 1197 

Then, in the strong confinement regime, the hypothesis of weak electric field limit should 
be valid as soon as the typical electric dipole interaction energy does not represent more 
than about 12% of the typical exciton confinement energy. If the micro-crystal radius 
is fixed at R = 10A, the highest electric field amplitude for which the weak field limit 
assumption stays acceptable is about _E™ ax pa 450kV.cm" 1 . Idem if the QD radius is fixed 
near R w 30A, the upper boundary of the strong confinement regime validity domain, 
the electric field amplitude inside the QD should not exceed E™ x m 16.7kV.cm~ 1 . These 
numerical results justify the choice of an electric field such that ps 12.5kV.cm _1 to 
compare theoretical predictions against experimental results, because it satisfies the 
weak field limit all along the strong confinement range of QD radius. 

In a more general manner, as soon as the semiconductor of the synthesized QD is 
chosen, the strong confinement regime domain of validity and the weak electric field limit 
condition consist of a set of two constraints, which should be optimized by choosing the 
QD radius and the electric field amplitude as functions of the Bohr radius, the Rydberg 
energy and the confinement energy of a trapped exciton. But, for future technological 
applications, this set of constraints will permit to determine conversely the best possible 
semiconductor for practical and technological reasons, by imposing the typical QD size 
and the order of magnitude of the maximal electric field amplitude to use. 



3.4-2. Case of computational data 
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Table 1. Stark shift for confined interactive electron-hole pair as a function of the QD 
radius including the Coulomb interaction and excluding the polarization energy, where 
terms scaling as oc cr 2 h are removed from the total Hamiltonian H mean value, i.e. 
C[ = and T^ r w —0.0333, in comparison with computational results [47] in spherical 
CdSo.12Seo.ss micro-crystallites. 



R 


A 


10 


20 


30 


40 


50 


A ^strong 

relative error 


meV 
mcV 


2.08 KT 4 
2.03 10~ 4 
<3% 


3.16 10" 3 
3.07 10~ 3 

<3% 


1.49 10" 2 
1.46 10^2 
«2% 


4.34 10" 2 
4.31 10~ 2 
<1% 


9.63 10" 2 
9.78 10~ 2 
<2% 



In the early nineties, a variational calculation on the same total Hamiltonian H 
with computational is performed, in order to study Stark effect in spherical micro- 
crystals [47] . The weak field limit is also considered and the Hamiltonian mean value is 
expanded in powers of the variational parameters <r ei h up to the second order. However, 
terms scaling as oc af h are neglected, while those scaling as oc a c a\ l are kept. Our 
approach suggests that both terms have the same order of magnitude and contribute 
to the electron- hole pair Stark shift. Leaving out such contributions implies that C[ 

should vanish in Eq. ( |T4l) . The expected Stark shift should be then affected, because 
3B'C" 

Tyar — C2 ~q — • Then, ps —0.0333 is independent of the semiconductor, and 

hardly represents about 20% of its value, when we account for contributions scaling as 

oc cr^ h , i.e. if C[ 7^ 0. Although the approximation made in [47] deeply changes the 

nature of Stark effect, and does not seem to describe correctly experimental results, 

except for very small QDs. Table [T] shows good agreement between our results if 

C[ = and computational ones from [47] , even outside the validity domain of the strong 

confinement regime, i.e. for R > 30A. This signifies that, in the strong confinement 

R 

regime, a first order expansion in — of the ground state energy Stark shift should be 

a* 

sufficient, at least when only the Coulomb interaction is included in the electron- hole 
pair Hamiltonian, because there is no particular constraint on the variational parameter 
a in [47]. 

4. Stark effects with polarization energy 

To investigate in more details Stark effects in semiconductor micro-crystallites and to 
especially integrate the electric dipole interaction point of view, we shall introduce the 
following polarization energy introduced in [8] to the total electron-hole Hamiltonian H 

1>1 K J 

where Pi(x) denotes the Legendre polynomial of index / and of variable x, and where 
we define, for a later purpose, the constants ati(e T ) as functions of the relative dielectric 
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constant e T as ai(e T ) = ^— — I EN*. We wish to apply a variational method 

K[le T + I + 1) 

to the new Hamiltonian H' = H + P, which takes into account the polarization energy 
P(r c , rjj. For this, we keep the form of the variational trial function $(r e , r^), since the 
polarization energy should not basically change the properties of the trial function. 



4-1. Strong confinement regime 

Here, the polarization energy mean value is expressed as an expansion in the variational 
parameters, which has the same form as the Coulomb potential mean value, where 
the constants A, B', C, C[ and C' 2 are replaced by functions of the relative dielectric 
constant e T . This is the reason why we adopt the same notations for these quantities, 
but with explicit dependence on e T , as shown in Table IA3I In fact, the polarization 
energy mean value in the quantum state defined by $(r e , r^) should be written as 

($|P(r c ,r h )|$> e 



2 



2 



— {A{e r ) + B'{e r )aR + C'{e { )a 2 R 

+C' 1 (s T )(a 2 e + a 2 h )R 2 + C' 2 (e v )a c a h R 2 + 0(a 3 P 3 )}. (18) 

In this formalism, we obtain expressions for the variational parameters and the Stark 
shift considering only the polarization energy or both the Coulomb interaction and the 
polarization energy in Stark effects by Eqs. (fT5l) and ( fTBT) in Section |3j by replacing all 
the appearing constants, in the first case, by the corresponding functions of e r and, in 
the second case, by the sum of both contributions. An important step in this calculation 
consists of the simple idea of rewriting the functions cti(e r ) in a such way that it becomes 
possible to perform analytically the summation of the series appearing in the expression 
of the polarization energy P(r e ,rh). As matter of fact, we remark that for / E N* 

Introducing A;(e r ) for / E N* as the relative error between the functions ai(e r ) and 
ai(e T ), the replacement of ai(e r ) by ai(e r ) is reasonable because it leads to negligible 
relative errors: e.g. for Cc/So.^Seo.ss micro-crystallites, for which e r = 4.0, we get that 
Ax(e r ) w 8%, A 2 ( £r ) w 3% and A,(e r ) = . £ *, 2l(1 \ u < 1%, for I > 3. The 

{s r + iy i{i + lj 

polarization energy is then written as 



e 2 e - 1 

P(r e ,r h ) 



f fogfl-4Vlogfl-4 



1-4 1 r h e r + i{ & V R 2 \ R 2 



2nRe r + 1 



(19) 



l>0 1>1 

where it is possible to sum a priori the series in the polarization energy because of the 
confining potential well. The mean value of P(r e ,rh) in the quantum state $(r e ,rh) 
contains only contributions of electron and hole such that r e ^ < R. 



Stark Effect of Interactive Electron-hole pairs in Spherical Semiconductor QDs 



19 



4-2. Comparison with experimental data 

The experimental parameters are those of Subsection 13.41 and the electric field amplitude 
inside the micro-crystal is kept at Ed = 12.5kV.cm~ 1 . 

4-2.1. Case of real experimental data 

We begin to evaluate the strong confinement regime validity region by considering Eq. 
(I17j) . on which we apply the appropriate changes according to the cases. If we only 
consider the polarization energy, the strong confinement regime is not actually valid, 
because it turns out that o~' < 0, in this case. Therefore, the interaction part of the 
variational function should be re i(r e h) = e' CT °' re . h . This reveals that the polarization 
energy is repulsive, while our approach is build on an attractive point of view for the 
interaction between the electron and the hole. If we consider both Coulomb potential 
and polarization energy, this problem no longer exists because the attractive effects of 
Coulomb potential are more important than the repulsive ones due to the polarization 
energy, so that a' remains positive. Thus, the strong confinement regime remains valid 
up to R < 1.0179a*, i.e. R < 50A for CdSo.12Seo.88 micro-crystallites. This domain 
of validity is twice as large as the one obtained if only the Coulomb potential is taken 
into account. In the same spirit, we evaluate an order of magnitude for the maximal 
electric field amplitude, for which the weak field limit should be assumed to be valid, 

as —7^ < 0.1186, when only polarization is taken into account, and —7^ < 0.1204, 

E c h E e h 

when polarization energy and Coulomb interaction are both included. From such 
numerical values, we should remark that the weak field limit is relatively insensitive 
to the interaction between the charge carriers. It may notably signify that this limit 
is relevantly chosen because it is an independent condition coming from the strong 
confinement regime validity domain. For example, this corresponds, for a QD radius of 
30A, to an electric field amplitude of about _E™ ax = 16.5kV.cm~ 1 , on one hand, and to 
an electric field amplitude of about £'™ ax = 16.8kV.cm~ 1 , on the other hand. 

As the electric field amplitude is fixed inside the CdSo.^Seo.ss m i cr( > cr y s tallites 
at Ed = 12.5kV.cm _1 , the weak electric field limit should be satisfied in the range of 
QD radii R < 30A and we can compare, at least in this domain, our predictions to 
experimental data. In the rest of the strong confinement regime validity domain, i.e. 
for radii 30A < R < 50A, the weak field limit is no longer valid. This may actually 
explain the significant divergence from experimental results in this region. Once again 
the analogy with rectangular quantum boxes should be quite helpful. In the opposite 
limit of strong electric field, the Stark shift undergone by the ground state of the electron- 
hole pair confined in a rectangular QD shows essentially a linear behavior in the electric 
field amplitude [40] , which is a totally different behavior in comparison to the quadratic 
dependence found here. Figure [2] clearly shows that the behavior of the polarization 
energy, if considered alone, is not satisfactory because it seems that its contribution 
does not counterbalance the divergence of the zeroth order Stark shift for QD radii near 
the upper boundary of the strong confinement validity region, while Figure [1] suggests 
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Figure 2. Stark shift for confined interactive electron- hole pair as a function of the 
QD radius including both the Coulomb interaction and the polarization energy up to 
the zeroth ( — ) or to the first order ( — ), where r°^ r sa —0.2045, and including only the 

polarization energy up to the first order ( ), where ss —0.0416, in comparison 

with experimental results (+) [34] in spherical C<i>S'o.i2>S'eo.88 micro-crystallites. 
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that the contribution to the Stark effects due to the Coulomb potential seems to consist 
of more important contributions. Figure [2] stresses this point of view. In the strong 
confinement regime, if we account for the Coulomb potential and the polarization energy, 
the results accurately fit the experimental data, except if the QD size begins to reach 
the lower boundary of the domain of this regime, in which the weak field limit is not 
valid anymore. The reason for this phenomenon is simple to understand. For such QD 
sizes, we remark that the first order term has the same order of magnitude as the zeroth 
order contribution to the Stark shift. Then, if it is reasonable to maintain a priori a 
perturbation point of view, the first order correction is not sufficient to describe correctly 
QCSE in the whole domain of strong confinement regime validity. It is perhaps advisable 
to continue the expansion to one or two further orders and to revise the definition of 
weak field limit. However, the computations become quite involved and we think that 
such approach does not really bring a significant improvement to the understanding of 
the Stark effects in QDs. 



4-2.2. Case of computational data 

Following the reasonings made in Subsubsection 3.4-2 . we add the polarization energy 
P(r c , r"h) to the total Hamiltonian of the electron-hole pair but impose that C[(e T ), C[ = 
0, in such a way that terms scaling as oc <j\ h do not contribute to the electron- hole ground 
state Stark shift. Once again, the new value for r^ r represents a small part from the 
one, in which all contributions are kept. Likewise, when we only consider the Coulomb 
interaction, Table [2] shows that there is still a good agreement between our results and 
computational ones from [47], over the whole validity domain of strong confinement 
regime. The divergence from experimental results is still significant after inclusion of 
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Table 2. Stark shift for confined interactive electron-hole pair as a function of the QD 
radius including the Coulomb interaction and the polarization energy where terms 
scaling as oc cr 2 > are suppressed from the total Hamiltonian H mean value, i.e. 
C|(e r ),Ci = and m —0.0446, in comparison with computational results [47] 
in spherical CdSo.12Seo.ss micro-crystallites. 
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A P 1 
i vr/ nuIO 

relative error 


meV 
mcV 


2.09 10" 4 
1.99 10~ 4 
«5% 


3.14 10" 3 
2.94 10" 3 
«6% 


1.45 10" 2 
1.36 10~ 2 
«6% 


4.11 10" 2 
3.90 10~ 2 
«5% 


8.93 10" 2 
8.55 10~ 2 
«4% 



polarization energy, even for small QD radius. This confirms that the terms scaling as 

oc (ip h play a relevant role in Stark effects and should not be discarded. While, this 

also legitimates the approximation ai(e T ) ~ on (e r ), because Table [2] suggests that this 

approximation does not seem over-estimate the polarization energy contribution to the 

Stark shift. Maybe, each term of the sum defining the polarization energy plays a role in 

Stark effect, but errors made term by term should not cumulate. A first order expansion 
R 

in — of the ground state energy Stark shift should not be sufficient and reinforces the 

a* 

idea according which it is necessary to carry on the expansion at least up to the second 
order. If the second order expansion does not improve the situation, this signifies that 
there should exist another reason for this divergence. Therefore, because the strong 
confinement regime validity domain is not affected be dropping of terms scaling oc <Tg h , 
it is reasonable to think that it comes from the failure of the weak field limit assumption 
for QD radius sufficiently close to its upper boundary. 

5. Conclusion 

By considering a simple EMA model under the assumptions of strong confinement regime 
by a infinite potential well and of weak electric field limit, we are able to obtain analytical 
results for Stark effect in semiconducting micro-crystallites with spherical shape. In 
the domain of validity of physical approximations, the numerical values we can deduce 
agree with experimental data. We furthermore clarify why other variational calculations 
predict numerical results, which markedly diverge from experimental ones. 

Despite these successes, our approach has been invalidated in a particular range 
of QD sizes, for which the strong confinement point of view still holds, but for which 
the weak electric field limit assumption fails. Thus, a future research work may focus 
on trying to apply the strong confinement regime in the limit of strong electric field 
indeed, even in a more general manner, to any electric field amplitude. The case of 
the weak confinement regime of the electron-hole pair is much more difficult, even in 
the weak field limit. Actually, in such a regime, the integration domain of integrals, 
we have to deal with, to compute the square of the trial function norm or the diagonal 
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matrix elements of physical operators consists of a half-rectangle instead of the domain 
V, which is a half-square. This implies an explicit break down of the electron-hole 
exchange symmetry. Then, a different approach may be needed. 

Appendix A. Constants 

In the following tables, we sum up all appearing constants and give their approximate 

f x dt 

values. The function Si(x) = / — sin(t) denotes the standard sine integral. 

Jo t 

Appendix A.l. Constants occurring in Stark effect expressions in absence of 
polarization 

Table lATI presents analytical expressions and approximate values of constants when only 
the Coulomb potential is taken into account. 

Appendix A. 2. Constants occurring in Stark effect expressions in presence of 
polarization 



We evaluate the polarization energy mean value using Eq. ( fl9i) . To this end, we compute 
integral representations of the polarization energy terms depending only on the radial 
coordinates r ej h of the electron and the hole following the reasoning, made in Subsection 
\3.2\ to get Eq. ([8]). This reasoning does not apply to the angular part. However, we are 
able to provide exact expressions for all the constants, which appears in the calculations, 
except for 5"', 7"' and 7"". For these, we obtain integral representations, which cannot 
be analytically computed at the moment. Their approximative values are computed 
numerically by using Wolfram Research Mathematica® 7. As exact expressions for 
other constants are quite cumbersome, we give only their approximate values. 

Let us define the constants (3' , /3", 7', 7", 7'", 7"", 5', 5" and 5"' by the expressions 

1-1% 1 r h 

= (3'- j'o-R + 5'a 2 R 2 + ^ (a c 2 + a 2 h )R 2 + 0(a 3 R 3 ), 

o 

-<*|log(l-f) + lo g (l-!)|*> 

= ,3" - YaR + S"a 2 R 2 + ^(ff 2 + a'QR 2 + 0(a 3 R 3 ), 

o 



Z>0 



2 _ Y'aR - 5"'a 2 R 2 + 2C<j a 2 + a A±^i \r> - ^LR 2 + 0(a 3 i^ 



2 < $ iEK^)' p ' (cos ^ hM 

i>i 
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Table Al. Definition, analytical expression and approximate value of constants when 
only the Coulomb interaction is taken into account. 



Name 



Expression 



Value 



Si(4?rJ 



.4 

Bx 
B 2 

D 
B' 

C 

a 



c 2 



D 2 
D 3 
D 
D' 
D" 
C" 



2 

25 



Si(27r) 



5 
8^ 



13 

24tt 2 



^Si(27r)- 



2 

3 
1 

Sl + f 
AB - 1 
1 1 



Si(47r) 



A(B 
B x 



2tt 2 
2AC 



37 



1207T 2 



12 
fh 
18 

2 13 147 

5 ~ 8^ + 64tt 4 
2 1 21 

15 ~ 8^ ~ 64tt 4 

320tt 4 2tt 5 \ 1 ; 2 
5£>i + lODg + D-3 
30 

3£>i + 4L» 2 + D 3 



5Z> 2 - -Da 



£>'■ 



45 

3D" - BC* 



1.7861 

0.6033 

0.2879 

0.6993 
0.2489 
0.2827 

0.0189 

-0.0339 

0.0160 

0.2589 

0.1173 

0.0710 

0.2539 
0.2195 
0.0115 
0.0187 



Table A2. Approximate value of constants appearing in ($|P(r c , rj 1 )|<i>). 



Name 


Value 


Name 


Value 


Name 


Value 


Name 


Value 


Name 


Value 




3.1144 


i 


2.3218 




1.3263 


6' 


0.9973 


5"' 


0.0533 


(3" 


0.7524 


7" 


0.5992 




-0.0704 


5" 


0.2708 







= _ Y" aR - 5 "'a 2 R 2 - -^R 2 + 0(a 3 R 3 ). 
1 18 V ; 

Table IA2I presents approximate values for constants which appear in the polarization 

energy diagonal matrix element (<&|P(r e , rh)|$), while Table R3l defines constants which 

appear in the polarization mean value Eq. ffl8|) and gives their approximate values in 

CdSo.uSeo.ss micro-crystals. 
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Table A3. Definition and approximate value of constants appearing in the polarization 
mean value Eq. (|18[) in CdSo,i2Seo.as micro-crystals. 



Name Expression Value for CdSo.i2>Seo. 



^ 2 — f- 2+ — 

1 F — 1 f F 



^Itif °' 5400 

«(e r ) i— l(j' + a'"-2C+-i— (i" + J'")l 0.4467 

A(e r ) ~/3(£r) -0.4467 

B'(e r ) - i 5(e r )S + 7 (e r ) -0.0993 

C'(e r ) -/3(e r )(B 2 -C)+ 7 (e r )B- ^ -0.0083 

C(( £r ) _^( er )-2^ (er )C 

C^( £r ) -^7?^ -0.0300 

18 
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